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Coulomb Gauge in One-Loop Quantum Cosmology 

The quantization program for gauge fields and gravitation in the presence of bound- 
aries is receiving careful consideration in the recent literature [1-12]. The motivations for 
this analysis come in part from the quantization of closed cosmologies [13-14], and in part 
from the need to understand the relation between different quantization techniques in field 
theory [15-16]. The choices to be made are as follows: (i) quantization technique; (ii) back- 
ground 4-geometry; (iii) boundary 3-geometry; (iv) boundary conditions respecting BRST 
invariance and local supersymmetry; (v) gauge condition; (vi) regularization technique. 
For a given choice of field theory, background and corresponding boundary 3-geometry, 
different quantization techniques and different gauge conditions have led to discrepancies 
for the semi-classical evaluation of quantum amplitudes [1-12]. This calculation has been 
performed within the framework of ^-function regularization, where the £(0) value yields 
both the scaling properties of the amplitudes and one-loop divergences of physical theories. 
On reducing a field theory with first-class constraints to its physical degrees of freedom 
before quantization [1-12], the resulting £(0) values disagree with the Schwinger-DeWitt 
A 2 coefficient [17-18]. This occurs both for compact Riemannian 4-manifolds without 
boundary [19-21] and for Riemannian backgrounds with boundary [1-12]. Moreover, fur- 
ther discrepancies have been found on studying the quantum theory of spin- 1/2 fields at 
one-loop about background 4-geometries with boundaries [5,22-25]. 

It therefore seems that the discrepancies found in the literature have at least two 
origins, as follows. 

(i) A geometrical source, owed to the singularity at the origin for manifolds with just one 
boundary. By this we mean that there is no regular vector field inside matching the normal 
at the boundary. Hence the normal and tangential components of physical fields inside are 
ill-defined, and it is impossible to achieve a consistent 3+1 split. Analogously, such a split 
cannot be obtained on compact 4-manifolds without boundary, if their Euler number does 
not vanish. 

(ii) A field-theoretical source, i.e. non-physical degrees of freedom and ghost modes yield 
contributions to £(0) which do not cancel each other on considering curved backgrounds 
and/or the presence of boundaries. 
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For example, for manifolds with just one boundary, the mode-by-mode analysis of 
Faddeev-Popov amplitudes, which relies on the expansion in harmonics of the potential 
(e.g. eqs. (4)- (5)) and hence on the eigenvalue equations for the modes multiplying each 
harmonic, disagrees with the Schwinger-DeWitt evaluation of the same amplitudes and is 
gauge-dependent. Moreover, on manifolds with two boundaries admitting a consistent 3+1 
split, the physical degrees of freedom of vacuum Maxwell theory, i.e. the transverse part 
of the potential, do not enable one to recover the full C(0) value [26]. These calculations, 
when combined with a previous analysis of discrepancies for spin-i fields [24-25], have 
motivated a mode-by-mode analysis of Euclidean Maxwell theory and linearized gravity 
[26-27] within the framework of Faddeev-Popov formalism [28-29] on manifolds with two 
boundaries. Remarkably, in this case the Schwinger-DeWitt and mode-by-mode formalisms 
for BRST-covariant Faddeev-Popov amplitudes are found to agree. While our previous 
analysis focused on relativistic gauges for such theories, it appears necessary to complete 
this investigation by considering non-relativistic gauges as well. Hence we here study the 
Coulomb gauge for vacuum Maxwell theory about flat Euclidean backgrounds with two 
boundaries, since such a gauge choice is more relevant for reduction to physical degrees of 
freedom and Hamiltonian formalism [30]. 

For this purpose, we use the version of ^-function technique [9,31] elaborated in refs. 
[6-8]. Following refs. [6-8], we write f n (M 2 ) for the function occurring in the equation 
obeyed by the eigenvalues by virtue of boundary conditions, and d{n) for the degeneracy 
of the eigenvalues parametrized by the integer n. One then defines the function 

oo 

I(M 2 ,s) = d ( n ) n ~ 2s log /n(M 2 ) . (1) 

n=no 

Such a function has an analytic continuation to the whole complex-s plane as a meromor- 
phic function, i.e. 

/(M 2 , s) = /pole(M2) + I R (M 2 ) + O(s) . (2) 
s 

The C(0) value is then obtained as [6-8] 

C(0) = /log + /pole(oo) - / P ole(0) , (3) 
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where 7i og is the coefficient of log M from I(M 2 ,s) as M — > oo, and I po \ e (M 2 ) is the 
residue at s = 0. 

To perform the one-loop analysis of vacuum Maxwell theory, we expand the normal 
and tangential components of the electromagnetic potential on a family of 3-spheres. With 
the notation of [2,9-10,26], one writes 

oo 

A (x,r) = J2Rn(r)Q {n) (x), (4) 

n=l 



A k (x,r) = J2 



n=2 



f n (r) SP(x)+g n (r)Pi n \x) 



(5) 



Of course, Q^(x), Sj^\x) and P^ n \x) are scalar, transverse and longitudinal harmonics 
on S 3 respectively [1-2,32]. 

Within the Faddeev-Popov quantization scheme [10,28-29], after adding to the Eu- 
clidean Lagrangian the Coulomb gauge-averaging term ^t -4 , one finds eigen- 
value equations on taking variations of the total Euclidean action with respect to the 
modes f n ,g n ,R n of eqs. (4)-(5) [10,26]. If a = 1, these equations take the following form 
(the decoupled mode will be treated separately): 

<J2 +-^--^ + Xn)fn(r) = 0, (6) 



where [33] 



dr 2 t dr r 

A n g n (r) + B n R n (r) = 0, (7) 
d n g n (T) + D n R n (r) = 0, (8) 



7 d 2 Id (n 2 - 1) A 

5 n = -(n 2 -l)(J- + i), (10) 
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(11) 




+ A 



■n 



(12) 



Moreover, in the Coulomb gauge, the ghost eigenvalue equation is found to be 




e n (r) = X n € n (r) Vra > 1 , 



(13) 



where the modes e n are the ones occurring in the expansion on a family of 3-spheres of the 
scalar field e(x,r) as Yl^Li e n( r ) Q^ n \x). Their contribution to £(0) is finally multiplied 
by -2 [9-10,26]. 

We here use magnetic boundary conditions [3,10,26]. Hence Ak(x, r) is set to zero at 



ingly the ghost modes e n [10,26]. 

The solution of eq. (6) is f n = A n I n (Mr) + B n K n (Mr), where M 2 = -X n . By virtue 
of our boundary conditions, such a solution should vanish at r = r_ and at r = r+, where 
T- , r + are the radii of the two concentric 3-sphere boundaries. Hence the corresponding 
eigenvalue condition can be written as 



By using the technique described in refs. [6-8] and outlined in eqs. (l)-(3), one can show 
that the contribution of transverse modes to £(0) is [26] 



From eq. (13) one can see that the only non-vanishing contribution from ghost modes 
is obtained when n = 1, which implies that A n = 0. The corresponding eigenfunctions are 
arbitrary functions of r which obey homogeneous Dirichlet conditions at r = r_ and at 
t = t+. The general eigenfunction can be written as 



the 3-sphere boundaries as well as the Coulomb gauge Q C {A) = r 2 A i , and correspond- 



I n (Mr-)K n (Mr+) - I n (Mr+)K n (Mr-) = . 



Cphys(O) 



1 



(14) 



2 




(15) 
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Note that we here face an entirely new situation, in that we have to deal with an infinite 
number of zero-modes. The inclusion of zero-modes into the general expression for ((0) 
was studied in mathematical papers [34] and for the problems of quantum gravity as well 
[35-36]. It is known that when we have a finite number of zero- modes we have simply 
to add it to the £(0) value, however, we do not know a priori what should one do with 
an infinite number of such modes. Hence we need the appropriate regular ization of an 
infinite number of our zero-modes. We try to achieve this being inspired by the ideas 
of (- regularization technique. For this purpose, we point out that the eigenfunctions 
(15) belong to the space whose elements can be parametrized by the natural numbers 
k = 1,2, .... All these eigenfunctions can be treated on equal footing. Thus, we are led 
to define the regularized dimension of this space as the regularized number of its basis 
elements, or, within the framework of ^-function technique (see below), as Cr(0) = — §, 
where (r(s) = J^^Li n ~ s * s the usual Riemann (- function, whose properties and values 
are well-known [37]. Thus 

Cghost(O) = -\ . (16) 

Of course, infinitely many definitions of regularized dimension are possible on considering 
the zeta-functions 

oo 
n=l 

where a is a real parameter. Our choice corresponds to the value a = 0, and it appears 
plausible to take into account that the sine functions of eq. (15) are labelled by the integers 
k = 1,2,3... only. 

By virtue of our boundary conditions, coupled gauge modes vanish at the boundaries, 
i.e. g n (j-) = g n (r + ) = 0, _R n (r_) = R n (r + ) = 0. The boundary conditions 

R n (r+) = R n (r_) =0 Vn > 1 

can be obtained from a non-trivial application of gauge invariance. In other words, one 
may start from a relativistic gauge condition written in the form (cf. [33]) 

$(A) = \A TrK + ^V l A l , (17) 
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where A is a dimensionless parameter and ( 3 )V l A; = t~ 2 A^ 1 . Now it is well-known that 
magnetic boundary conditions in relativistic gauges imply that any such should 
vanish at the boundaries [26]. When combined with homogeneous Dirichlet conditions on 
the tangential components of the electromagnetic potential, which require that 

f n (r+) = / n (r_) = , g n (r+) = g n {r_) = , Vn > 2 , 

this implies that Aq(x,t + ) = A (x,T-) = 0, i.e. 

R n (r+) = R n (r-) =0 Vn > 1 . 

If gauge invariance is respected in quantum theory, this set of conditions holds for all 
values of A. Hence, on taking the limit as A in (17), one recovers the Coulomb gauge 
we are interested in, subject to the same boundary conditions on g-modes and .R-modes, 
providing gauge invariance holds. The legitimacy of this procedure depends crucially on a 
direct proof of gauge invariance in the presence of boundaries, not relying on the formal 
arguments frequently presented in the literature, and is the object of a paper in preparation 
by ourselves and other co-authors [33] . We are then able to obtain, in particular, boundary 
conditions on the decoupled mode which would otherwise remain totally arbitrary. 
Moreover, eqs. (7)-(12) imply that 



Rn(r) = 



9n{r) 



[(n 2 - 1) +M 2 r 2 



(18) 



where g n {r) obeys the second-order differential equation 



(n 2 - 1) + M 2 r 2 
(n 2 - 1) +M 2 t 2 



M 2 T 2 d 2 g n 



dr 2 
2 (n 2 - 1) 



3(n 2 - 1) +M 2 t 2 



-9n 



M 2 T d9n 

dr 



(n 2 - 1) + M 2 t 2 



9n = . 



(19) 



Note that, if M = 0, the limiting form of (19) is 



(n 2 - 1) 



g n = Vn > 2 , 
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whose only solution is g n (r) = 0, Vn > 2 and Vr G [r_, r+]. If M does not vanish, we may 
regard (19) as a second-order ordinary differential equation. Remarkably, one then deals 
with an overdetermined problem, in that both g n and g n have to vanish at the boundaries 
in the light of boundary conditions and of eq. (18). Hence the only solution is the trivial 
one, i.e. g n (r) = Vr <E [t_,t+] and \/n > 2, and similarly for R n (r),\/n > 2. Thus, 
coupled gauge modes give a vanishing contribution to the full C(0) : 

Ccoupled(O) = . (20) 

By contrast, the decoupled mode Ri(t) can be represented by an arbitrary function 
vanishing at r = r_ and at r = r+. Hence, one can apply again the argument leading to 
eq. (16), which implies 

Cflx(0) = ~. (21) 
By virtue of eqs. (14), (16), (20)-(21) one finds 

C(0) = Cphys(O) + Ccoupled(O) + ( Rl (0) - 2C ghos t(0) = . (22) 

One can also consider another approach to boundary conditions for the electromag- 
netic field subject to the Coulomb gauge. For this purpose, we can view the normal 
component A of the electromagnetic field as a Lagrange multiplier which should not be 
included in the gauge conditions and which should be integrated over not only inside the 
manifold under consideration, but also on its boundaries [30]. In this case the homoge- 
neous mode R\ should be excluded, since it does not correspond to any constraint (really, 
the constraint ViF l ° of the electromagnetic field has no homogeneous modes). Moreover, 
homogeneous ghost modes should also be discarded, because no homogeneous modes cor- 
respond to the Coulomb gauge condition (the problem of discarding ghost zero-modes was 
discussed in a slightly different framework in [38]). 

Thus, in such an approach to the problem of boundary conditions and the treatment 
of Lagrange multipliers and ghost modes, there are no non-zero contributions to the full 
C(0) from the decoupled mode and from ghost zero-modes (cf. eqs. (16), (21)). However, 
in this case we have a non-trivial contribution from the coupled gauge modes. Hence 

8 



Coulomb Gauge in One-Loop Quantum Cosmology 



our second-order equation (19) is no longer over determined, since there are no boundary 
conditions on R n (see (18)), and we have a Dirichlet boundary- value problem involving 
(7 n -modes only. One can then show that the basis function of eq. (19) can be represented 
as 

9n(r) = w v {Mt) v n {Mr) x n (M, r) , (23) 

where w v is a linear combination of modified Bessel functions I v and K v with v 2 = 2(n 2 — 1). 
Then 



Vn(Mr) = 
and x n (M, r) obeys the equation 



[(n 2 - 1)+M 2 t 2 



M 2 T 2 



d 2 x 



" + - ^ + F(n, M, r)x r 



dr 2 t dr 



+ 2 



dw v dx n 
dr dr 



(24) 



where 



F(n, M, t) 



(n 2 ~ 1) 
M 2 r 4 



(^1) \ ^ / 3(n»-l) + MV 
1 ^ \ ^ M 2 r 2 / V (n 2 - 1) + M 2 t 2 



(n 2 -l) ^ + (n 2 -l)' 



M 2 



V ^ M 2 r 2 7 



x 



x 



3 (n 

7 + V 



M 2 r 6 



1)/ fn 



M 2 r 2 



(25) 



The solution of (24) has the following asymptotic forms: 



x n (M, t) ~ A n + as M -»• oo , 



where A n is a constant, and 



x n (M, r) ~ exp 



(n 2 - 1) 
Mr 



as M -»• 



Hence one can show that the only function in the product (23) which can give a non- 
trivial contribution to £(0) is u^, and its contribution can be easily calculated along the 
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lines described in refs. [6-8,24-27] and applied above to the calculation of the contribution 
of physical modes (eq. (14)). This leads to 

Ccoupled(O) = \ . (26) 

Combining eqs. (14) and (26) one has again 

C(0)=Cphys(0)+Ccoupled(0)=0. 

Thus, our result obtained with another treatment of boundary conditions coincides with 
the one resulting from eq. (22). 

Interestingly, our result for the full ((0) agrees with previous calculations for relativis- 
ts gauges [26] and with the A 2 coefficients of the Schwinger-DeWitt technique [26]. For 
our manifold with two 3-sphere boundaries, the general formulae for the Schwinger-DeWitt 
A 2 coefficient [3,5,18] yield zero, since the contributions of the two boundaries cancel each 
other. 

We should emphasize that part of our analysis relies on the Riemann-C regularization 
of the infinite dimension of the space of zero-modes for the decoupled component of A and 
for the ghost field. We can only say that this procedure leads to the remarkable cancellation 
occurring in (22), but we cannot as yet put it on more rigorous grounds. Moreover, our 
boundary conditions in the Coulomb gauge result, in the first approach, from a non-trivial 
limiting procedure out of the well-established set of boundary conditions imposed on study- 
ing relativistic gauges. Such a limiting procedure leads to a non-vanishing contribution 
to £(0) resulting from ghost modes. As it is well-known, this non-vanishing contribution 
is absent in the case of flat space without boundaries. The ultimate justification of our 
approach can be obtained by proving in a direct way the gauge invariance of quantum 
amplitudes for problems with boundaries, rather than assuming it. A detailed analysis of 
such a problem can be found in ref. [33]. It seems encouraging, however, that boundary 
conditions for the Coulomb gauge derived from the familiar, non-relativistic framework, 
lead to a C(0) value equal to the one in (22). 
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In agreement with recent work by the authors [26-27], our analysis of the Coulomb 
gauge adds evidence in favour of gauge- and ghost-modes being necessary to obtain gauge- 
invariant amplitudes for manifolds with boundaries. In general, there is no exact cancel- 
lation of such contributions on non-trivial backgrounds. Moreover, we should say that 
a very recent paper by Moss and Poletti [39], relying on previous work by Vassilevich 
[40], has corrected all previous Schwinger-DeWitt calculations for spin> for manifolds 
with boundaries. Hence they find full agreement between our investigations of the con- 
formal anomalies for spin-i and spin-1 fields [22-26], and their corrected analysis, even 
in the 1-boundary case. However, the problem remains of understanding why, in the 
1-boundary case only, the mode-by-mode evaluation of one-loop quantum amplitudes is 
gauge-dependent, as shown in [26]. As far as we can see, this last remaining discrep- 
ancy seems to point out to some serious limitations which apply when the background 
4-geometry does not admit a well-defined 3+1 decomposition. 

The problem of the correct choice of technique for one-loop calculations in quantum 
cosmology is also important when one investigates the properties of the wave function 
of the universe [13-14]. These properties are crucial for studying the relation between 
quantum gravity, inflationary cosmology and particle physics [41-42]. 
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